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In an attempt to look for a viable mechanism leading to a present day accelerated expansion, we 
investigate the possibility that the observed cosmic speed up may be recovered in the framework of 
the Rastall's theory, relying on the non - conservativity of the stress - energy tensor, i.e. T^.^j 7^ 0. 
We derive the modified Friedmann equations and show that they correspond to Cardassian-like 
equations. We also show that, under suitable assumptions on the equation of state of the matter term 
sourcing the gravitational field, it is indeed possible to get an accelerated expansion, in agreement 
, with the Hubble diagram of both Type la Supernovae (SNela) and Gamma Ray Bursts (GRBs). 

v,^ ■ Unfortunately, to achieve such a result one has to postulate a matter density parameter much larger 

than the typical — 0.3 value inferred from cluster gas mass fraction data. As a further issue, 
we then discuss the possibility to retrieve the Rastall's theory from a Palatini variational principle 
I ■ approach to f{R) gravity. However, such an attempt turns out to be unsuccessful. 

o 



I. INTRODUCTION 



The observed cosmic speed up [l|, H, 0] questions the validity of General Relativity (GR) on large scales. In fact, 
if on one hand the model of gravitational interaction as described by Einstein's theory is in agreement with many 
observational tests on relatively small scales, as Solar System and binary pulsars observations show it is well 
known that in order to make GR agree with the observed acceleration of the Universe the existence of dark energy, 
a cosmic fluid having exotic properties, has been postulated. Actually, many candidates for explaining the nature of 
dark energy have been proposed (see e.g. Q, and references therein), some of them relying on the modification 
of the geometrical structure of the theory, some others on the introduction of physically (up to day) unknown fluids 
\ into the equations governing the behaviour of our universe. In this context, we want to consider here a generalization 
of Einstein's theory based on the requirement that the stress-energy tensor for the matter/energy content is not 
conserved, i.e. Tj^.^ ^ 0, to check if it is able to describe the accelerated expansion of the universe. In particular, we 
focus on the so called Rastall's model (7j], and we investigate the conditions that the parameters of the theory have to 
\^ fulfill in order to reproduce the observed accelerated expansion of the Universe. The field equations of this model were 
5^ , deduced by Rastall as ad hoc modifications of GR ones, so as to have a non conserved stress-energy tensor; however, 
there have been subsequent attempts to deduce Rastall's field equations from a variational principle, but none of them 
have succeeded [1, Since these field equations are similar to those that in f{R) gravity are obtained by a Palatini 
variational principle from a Lagrangian depending on an arbitrary function / of the scalar curvature R (lol . , we 
aim at verifying whether Rastall's model can be considered a particular case of .f(R) grav ity. Furthermore, prompted 
by a recent paper |13] , we check the possibility that Cardassian-like equations [ij, [lj| can be really derived from 
Rastall - like theories. 

The plan of the paper is as follows: we will firstly give an introduction to Rastall's model in Sect. 2, while the 
corresponding cosmological scenario and the analogy with the Cardassians expansion model is worked out in Sect. 3. 
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In order to check the possible viability of the Rastall's proposal, we fit the model to the SNela and GRBs Hubble 
diagram as detailed in Sect. 4, while Sect. 5 focuses on a more theoretical aspect investigating the possible connection 
with f{R) theories in the Palatini formalism. Conclusions are finally presented in Sect. 5. 



II. RASTALL'S MODEL 

In 1972 P. Rastall 0] explored a model in which the stress - energy tensor of the source of the gravitational field, 
T^i^, was not conserved, i.e. the condition T^^ 7^ is imposed a priori. 
Indeed, Einstein equations^ read 



Gf_ii, — Rf_ii, — -^g^vR — i^grT^iv , (1) 

where the Ricci tensor is obtained from a metric connection, so that i?^,^ — Rf^^, (g) and the scalar curvature R has to 
be intended as i? = R{g) ~ g^^Rapig); furthermore we have set kgr = 

These equations naturally imply the stress - energy tensor conservation as a consequence of the contracted Bianchi 
identities, 

= . (2) 

It is therefore worth wondering whether it is possible to fulfill the requirement Tjf.^ ^ without violating Eqs.©. 
A possible way out could be introducing further geometrical terms on the right hand side of Einstein equations, but 
one should ask whether this makes sense. Actually, if we insist in deriving these relations from a metric variational 
approach, the sudden answer would be of course negative: in this case the stress - energy tensor would be surely 
conserved by construction, so no way to escape the conditions T/^^ = 0. 

Another remark against the non - conservativity focuses on the equivalence principle: as a matter of fact, the 
conservation of the stress - energy tensor is tested with high accuracy in the realm of Special Relativity (SR). However, 
one has to go easy with such an approach, as this principle could be misleading [l5| . To give an example, when passing 
from GR to SR, we completely miss the information provided by terms explicitly depending on the curvature tensor, 
R^g^, as it becomes identically zero in Minkowski spacetime. This means that the two sets of equations 



(3) 



and 



V.r + i?^;3.V„j''VV = 0, (4) 

give exactly the same equations, i.e. daj" — 0, in SR. So, the straightforward application of the equivalence principle 
in writing conservation laws should be carefully considered. 

The question is now how to pick up a proper geometrical term such that the Bianchi identities are still not violated, 
but nevertheless the conservation of the stress - energy tensor of the gravity source is violated. To resume, we ask for 
a four-vector, say a^, such that (i) Tlf.^ = a^] (ii) ^ on curved spacetime, but = on flat spacetime in order 
not to conflict with the validity of SR. Both these properties hold for the Rastall's proposal, that is 



Ti:,^ = Ai?,, , (5) 

A being a suitable non -null dimensional constant. 

Because of the assumption ([5]) , the field equations are obviously modified and now read 



(6) 



^ Throughout the paper, spacetime is assumed to have the signature {-|-, — , — , — ), and Greek indices run from to 3. 
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where is a dimensional constant to be determined in order to give the right Poisson equation in the static weak - field 
limit. It is manifest that in vacuum, where = 0, Rastall's field equations © are equivalent to GR ones. 

As a matter of fact, the same set of equations can be obtained as the result of guesswork, that is assuming the left 
hand side of the sought after equations to be a symmetric tensor consisting only of terms that are linear in the second 
derivative and/or quadratic in the first derivatives of the metric 16]. Moreover, the time-time component of such 
equations must give the Poisson equations back for a stationary weak- field. Accordingly, the only requirement we 
drop with respect to the derivation of the Einstein equations is the one concerning the conservation of T^^. Hence, 
starting from G^^, = CiR^^, -\- C2Rg^vi with Ci and C2 appropriate constants, we end up with Eqs.® again provided 
that we set : 



2(3- 
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-2Ci) ' 



(8) 

where we have chosen to rewrite all of them in terms of the constant Ci . Note, in particular, that the coupling constant 
between matter and geometry, k^, is not the same as in GR, unless Ci = 1, that is A = (i.e., we consistently go 
back to GR). 

Taking the trace of Eqs.® gives us the structural or master equation : 

{AKrX - I) R ^ KrT . (10) 

For a traceless stress - energy tensor, T = (as for the electromagnetic tensor) and two possibilities arise. The first is 
that i? = so that we get no difference with standard GR. On the other hand, one could also solve Eq. p^ setting 
Kr\ =1/4, whatever the value of R is. However, inserting this condition in Eq.®, we get a complex value for Ci 
which is clearly meaningless. Therefore, we hereafter assume that k^A ^ 1/4. 

Another fundamental question concerns geodesic motions. As it is well known, the equations Tj^.^ — are nothing 
but the equations of motion of the fluid we are dealing with. The problem is then, what sort of curves are described 
in a curved spacetime by a fluid whose stress - energy tensor is not conserved ? Following the calculations made by 
Rastall, we find that in his model geodesies are those curves characterized by the fact that the scalar curvature R is 
constant along them. Moreover, it is still possible to speak of conservation of energy for an ideal fluid [7|, but again 
provided that R is constant along the time -like four - velocity vector of the fluid, it^. The question remains whether 
particles creation takes place in the regions where this condition does not hold. 

It is worth mentioning that the Rastall's equations ^ can be recast into the same form as the usual Einstein ones. 
Indeed, one can immediately write 



G — KrSfiu , (11) 

where 



By construction, this new stress - energy tensor is conserved, S^.^ = 0. On introducing 5'^;^, we can recover all the 
known solutions of Einstein GR by simply taking care of the difference between 5*^^ and T^jy. Furthermore, if we 
assume T^j^ — {p + p)u^Uy — i.e. the source is a perfect fluid with energy density p and pressure p, we can 

explicitely work out an expression for This turns out to be still a perfect fluid provided we redefine its energy 
density and pressure as 



PS = 

PS = 



(3Kr A — l)p + SKrXp 

4:KrX - 1 
KrXp + (KrA — l)p 
4:KrX - 1 



(13) 
(14) 
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In order to obtain the value of the couphng constant Kr , we remember that the time - time component of the modified 
equations should recover the Poisson equation in the static weak -field limit. One thus gets : 

3Kr\ - ) = KGR , (15) 



4k^A- 1 V 2^ 

whence it is immediate to derive exactly the same coupling as in Einstein gravity only when A = 0, that is when the 
conservation of T^^ is granted. 

III. A CARDASSIAN ANALOG AND THE COSMIC SPEED UP 

It has been recently claimed that a Cardassian- like [H, [ij modification of the Friedman equation in the form 

H' = ^[p + Bit){p-3pr] , (16) 

can be obtained from Rastall-like equations, where B{t) is a function of the cosmic time t. We would now like to 
show that, although it is indeed possible to recast the Rastall's theory equations in such a way that a Cardassian- like 
model is recovered, the parameter B in must be a constant. 

To this aim, we derive the cosmological equations for the Rastall's theory. We first remember that, when the 
isotropic and homogenous Robertson - Walker (RW) metric is adopted, in GR one gets the usual Friedmann equations 

- , (17) 

H + H' = -H^ip + Sp) , (18) 
6 

where H = a/ a is the Hubble parameter, a the scale factor and a dot denotes the derivative with respect to cosmic 
time t. To get the corresponding equations for the Rastall's theory, one has to insert the RW metric into Eqs.® and 
consider the only independent equations that can be obtained, that is 

3a^ — 6KrA (a^ + aaj — KrO^ p , (19) 
+ 2aa + Gk^A (a^ + aa) ~ —UrO^p , (20) 

respectively. The master equation thus becomes 

6 (a^ -I- aa) = —Kra^ {p — 3p) , (21) 
so that multiplying Eq. ([20|) by —3 and then adding to Eq. (fT9|) we finally get the first modified Friedmann's Equation: 

a _ Kr p + 3p-6KrX{p + p) ^^2) 
a 6 4:KrX — 1 ' 

which makes it possible to directly infer the sign of the acceleration. To obtain the second modified Friedmann's 
equation, it is easier to proceed in a slightly different way. Let us first take the Rastall's equations in the form 

i?,. = -.(r,.-i|^rg,.) . (23) 
By inserting the RW metric and adding up Ea. (jl9p with three times Ea. (j20p . we eventually obtain 
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2KrX — 1 , 

ip + 3p) + —^^ip^3p) 



(24) 
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It is then only a matter of algebra to rearrange Eqs. ((22l) and ((24l) to write them as 



-iKrX - 1 



(25) 



ij^^^^^;f^b + .-4..A(p-p)] , (26) 

which reduce to the standard Friedmann equations P?)) and when the parameter A is switched off. Note also that 
Eq. ([25ll has indeed the same expression as the Cardassian- like Eq.(fT6|l provided we set n = 1 and accordingly redefine 
the parameter B. However, it is straightforward to show that B must be a constant. Indeed, from the condition 

S'' = , (27) 



it is immediate to demonstrate that Hj-X must be a constant by simply inserting the master equation ((2T|) into Eq. ((27l) 
and using the Rastall's requirement = XR^. So, an equation like (|16p . cannot be self-consistently obtained in 
Rastall's model. 

Moreover, since T^^ is a perfect fluid and remembering the definition of 5*^^, Eq. (|27p . we get 

p + 3Hip + p)^ (p-3p) , (28) 

which generalizes the continuity equation for the Rastall's theory. 

It is worth noticing that, even without integrating the equations, one can immediately predict whether the universe 
expansion is accelerating or not by simply studying the sign of the right hand side of Eq. (j22p . 

Assuming for simplicity that the equation of state of the perfect fluid is a constant, i.e. settings — wp, the condition 
a > selects two possible regimes, the first one being 

^'^rX - 1 , , 

" > 3(1 - 2«.A) ' 

provided k^A > 1/4. When A = 0, however, the above relation reduces to w > —1/3 in contrast with the GR result. 
We have therefore to choose the other solution, namely 

6KrX-l 

"^ 3(1-2^.A) ' 

with KrX < 1/4. The right hand side of (|30p may be positive or negative depending on the value of k^A. More 
precisely, if 1/6 < KrX < 1/4, then the rhs is positive, while it is negative for < KrX < 1/6. It is worth stressing 
that, however, the model always gives a monotonic behaviour: always decelerated or, as in the above analyzed case, 
always accelerated. 

By the way, in the spirit of Cardassians, the only sources of gravity are radiation and matter. In particular, the 
recent epoch is driven by the matter content, described as a perfect fluid with equation of state w — 0. With this 
constraint, it is easy to show that an accelerated behaviour is obtained for i < KrX < j, whereas we have a decelerated 
expansion choosing the following values KrX < ^ or k^A > |. 

IV. RASTALL'S MODEL CONFRONTED WITH THE DATA 

Neglecting the radiation component, the only fluid sourcing the gravitational field is the standard matter, which 
can be modeled as dust, i.e. p = 0. In such a case, the continuity equation (|28p is straightforwardly integrated giving : 

p cx po(l + z)'''^" , (31) 
with a subscript denoting present day quantities, z — 1 = i (having set oq — 1) the redshift, and 
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an effective equation of state (EoS) for the dust matter, from which a relation between Weff and k^A is easily deduced. 
Note that, for A = 0, one recovers the usual matter scaling p cx (1 + z)^, while deviation from the standard behaviour 
occurs in A 7^ Rastall's theory. Such a different scaling is not surprising at all being an expected consequence of the 
non - conservativity of the stress - energy tensor. Inserting back Eq. (PT|) into Eg.ipS)). we get 

= H'^/Hl = (1 + zf""' , (33) 
which is all what we need to compute the luminosity distance 



1 

DL{z,Weff,h)=dH{l + z) ^^dz\ (34) 

with the Hubble radius dn ~ c/Hq ~ 2>h~^ Gpc and h the Hubble constant in units of 100 km/sMpc. We have now 
all the main ingredients to test the viability of the Rastall's model by fitting the predicted luminosity distance to 
the data on the combined Hubble diagram of SNela and GRBs. To this aim, we maximize the following likelihood 
function : 



C{weff, h) (X exp 



2 2 
XsNela + XgRB 



X exp 



<^HST 



(35) 



with 



J^SNela 

XsNela = 



^J■obsiz,) - Hth{Zi,Weff,h)' 



(36) 



Xgrb 



J^GRB 

E 

i=l 



Hobs (Zi) - IJ,th{Zt, Weff , h) 



(37) 



The terms in (|35p take care of the Hubble diagram of SNela and GRBs, respectively, and rely on the distance 
modulus defined as 



IJ.th{z,Weff,h) = 25 + 5\og DL{z,Weff,h) 



(38) 
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We use the Union [l8| dataset for SNela and the GRBs sample assembled in Cardone et al. |l9| to set the observed 
quantities {fJ-obs,o'i) for the SNela and GRBs, respectively. Since the Hubble constant h is degenerate with the 
(unconstrained) absolute magnitude of a SN, we have added a Gaussian prior on h using the results from the HST 
Key Project ^20] thus setting {hHST,crHST) = (0.72,0.08). 
The best fit model turns out to be 

(zx;e//,/i) = (0.55,0.68) , 

giving 

xlNeiJd.o.f. = 1.08 , xlRB/d.o.f. = 2.07 , 

where d.o.f. = NsNeia + Nqrb ~ Np is the number of degree of freedom of the model, with NsNeia = 307 the 
number of SNela in the Union sample, Ngrb = 69 the number of GRBs, and Np — 2 the number of parameters of 
the Rastall's model. While for SNela we get a very good reduced x^, this is not the case for GRBs so that one could 
be tempted to deem as unsuccessfuU the fit. Actually, Fig. [1] shows that the model is indeed fitting quite well both 
the SNela and GRB data so that the large value of Xq^^/c?.o./. should be imputed to the large scatter of the high 
redshift data around the best fit line, not taken into account by the statistical error on the GRBs distance modulus. 
In order to further test the model, one can consider the constraints on the matter density parameter. Since our model 
only contains matter, one could naively think that — 1- Actually, one must also take into account that VIm is 
defined using the GR coupling constant kgr which is related to the Rastall's coupuling Kr through Eq. p5)) . It is then 
a matter of algebra to show that = Wefj so that, after marginalizing over h, we get the following constraints : 

nM = 0.55^°02 +"-05 



X 



'M — U.OO_Q Q3 _o.05 ' 

where we have used the notation x'^iy^^y^ to mean that x is the median value of the parameter and {x + xi, 
yi),{x + X2,x — 2/2) are the 68% and 95% confidence ranges respectively. Note that the value thus obtained is in 
strong disagreement with the typical Om — 0.3 obtained from both the cosmic microwave background radiation data 
and cluster gas mass fractions. Such a large matter density parameter is clearly unacceptable and represents a strong 
evidence against the Rastall's model. It is worth noting that such a result could be qualitatively foreseen considering 
that, because of the non- conservativity of the matter stress - energy tensor, a sort of matter creation takes place thus 
increasing fi^/ and leading to the final disagreement. 

V. RASTALL'S FIELD EQUATIONS FROM A PALATINI APPROACH 

Rastall's model has been initially motivated by the need for a theory able to allow a non conservativity of the 
source stress - energy tensor without violating the Bianchi identities. As such, the original theory was based on purely 
phenomenological motivations starting directly from the field equations without any attempt to derive them from 
a variational principle. We intend here to readdress this issue looking for a connection to f{R) theories of gravity, 
which are obtained by a variational approach, starting from a gravitational Lagrangian depending on a function / of 
the scalar curvature R (see [13, [ill references therein). These theories can be studied in the metric formalism, 
where the action is varied with respect to metric tensor, and in the Palatini formalism, where the action is varied 
with respect to the metric and the affine connection, which are supposed to be independent from one another. In 
particular, we want to verify here if Rastall's theory can be considered a particular case of Palatini f{R) gravity, 
which would suggest how to obtain it from a variational principle. 

To this aim, let us first rewrite the field equations in Rastall's theory as follows : 



-Rpi. - T^RgtLv = i^r {Tp„ - >^R9fiiy) , (39) 



where we stress that i?^^ = Rf^^{g) and R = R{g) — q"^ RaB(ci) for reasons that will be clear soon. Assuming that 
f'{R) — df{R)/dR ^ 0, the corresponding field equations [l3;[2l| for a /(i?) theory in the Palatini approach can be 



written as 



^-m9) 



9^.u , (40) 
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(41) 



where i?^i/(r) is the Ricci tensor of a torsionless connection F independent of the metric g, R{T, g) — g'^'^Rf_ii, is the 
scalar curvature buih with the metric tensor and we have denoted with R(^u) the symmetric part of R^v 

To compare Eqs. (p9)) and (|40l) - (I4ip . we first rewrite Ri_iv(T) in terms of the Ricci tensor of the metric, Rfj_^{g), 
which gives [H, [^J] 



(42) 



where we have dropped the various dependences for sake of clarity. Note that we have here introduced the vector field 



which allows to rewrite Eqs. ((4T|) as 



-V„ [ln/'(i?)] 



(43) 



and to define the rank two tensor field 



(44) 



(45) 



By saturation of indices in Eqs. (l42p . one naturally obtains the corrensponding relations between the scalar curvature 
R{T,g) of the affine connection and that of the metric, R{g), reading: 



R^R + 2V„6" + V"6„ + . 

2 



(46) 



Then, inverting Eqs. (|42|) and ([46]) and introducing into the Rastall's equations p9|) the Ricci tensor and scalar 
curvature of the metric in terms of the corresponding barred quantities, we obtain the following relations : 



R^u ^g^uR — ^rT^u ~t~ l^r^Qiiu 



(47) 



Note that the structural equation from (HT)) reads 



R^ , : r + 2V„6" + V"6„ + -b„b" . (48) 

Considering Eq. ((46l) , it is immediate to see that Eq. (l48l) is the same as Eq. ((TO)) thus showing us that the same relation 
between the metric scalar curvature and the trace of the source stress - energy tensor holds both for the /(i?) and 
Rastall's theories. 

From the very definition of ba and a manipulation of Eq. (|^D|) . it is easy to show that, for a general f{R) theory, 
the following result is valid 



__ f"{R) 
" - ""GR _ . 

This means that in empty space, or in case of a traceless stress - energy tensor, the vector field 6q is null and, according 
to (j4T|) . the connection F is metric, so that R^,^{V) = Rf^,, and R{T,g) — R. As a consequence, the field equations 
become 
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/'(i?) '^"^ 2 



f{R) 



R 



(50) 



In such a case, we can easily equate the rhs of Eg. ([50)1 with that of Eg. ([55)1 and then work out an analytical expression 
for f{R) such that the two sets of field equations are identical. Provided k^A 7^ 1/2 (since it would correspond to 
f{R) = 0) and setting a = 2Kr\ — 1, we finally get 



f{R) = foR-i . (51) 

Note that A = gives a = —1 hence f{R) = R and we correctly recover the GR f{R). On the contrary, for A 7^ 0, 
inserting Eg. ([5T|l into Eg. (|Sn|) makes it possible to write the field equations as 

Rfj.i' ^ 2^-91^1' — yjjj^'Tf^i' ~ KrXRg^i, . (52) 

Comparing Egs. ([55)) and shows that, although having the same structure, they are different because of the diverse 
couphng, being kgr 7^ Kr- Remembering Eq.®, we see that the only way to make Eos. ([55|) and to coincide is 
setting 



f'iR) = C, (53) 

which is clearly incompatible with the solution ((5T|) . unless i? is a constant. We therefore conclude that the Rastall's 
theory may be recovered from a f{R) gravity model in the Palatini formalism only when it reduces to GR. Put in 
other words, this means that Palatini f{R) gravity can not provide a viable Lagrangian from which the Rastall's field 
equations may be derived. 



VI. CONCLUSIONS 



In this paper we have focused on Rastall's theory of gravity, whose field equations are based on the request of 
the non - conservativity of the stress - energy tensor, i.e. T,^^ 7^ 0. In particular, we have reexamined this model of 
gravity to investigate the possibility that it could reproduce the observed cosmic speed up. To this end, we have 
explicitly worked out the modified Friedmann equations and confronted the model predictions with the available data 
concerning type la Supernovae (SNela) and Gamma Ray Bursts (GRBs) : what we have showed is that it is possible 
to get an accelerated expansion that is in agreement with the Hubble diagram of both SNela and GRBs, even if there 
is unfortunately no possibility to reproduce an accelerated-decelerated-accelerated expansion for our universe as it 
seems to be requested. These results have also a major drawback: indeed, to get them it is necessary to postulate a 
matter density parameter much larger than the typical ^Im — 0.3 value inferred from cluster gas mass fraction data. 
As a consequence, Rastall's theory is not in agreement with current cosmological observations. We have also shown 
that Cardassian-like modifications of Friedmann equations are obtained in Rastall's model but, contrary to recent 
claims, they cannot contain time-dependent parameters. Finally, in an attempt to derive Rastall's field equations 
from a variational principle, we have studied the connection with f{R) gravity in the Palatini formalism, and we have 
concluded that the Rastall's theory may be recovered from a f{R) gravity model in the Palatini formalism only when 
it reduces to GR; thus f{R) gravity cannot provide a variational principle for deriving Rastall's field equations. 
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